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Planar  solid  oxide  fuel  cells  with  yttria-stabilized  zirconia  electrolytes  typically  operate  at  tempera¬ 
tures  in  the  range  of  700-850  °C.  The  maximum  temperature  is  limited  by  the  use  of  ferritic  stainless 
steel  interconnects  which  offer  significant  advantages  such  as  low  cost  and  high  thermal  and  electronic 
conductivity  over  traditional  ceramic  interconnects.  However,  these  alloys  rely  on  the  formation  of  a 
protective  chromia  scale  for  oxidation  resistance  that  results  in  an  increase  in  ohmic  resistance  and  can 
volatilize  leading  to  a  loss  of  cathode  catalytic  activity.  To  better  understand  the  oxidation  behavior  of 
chromia  forming  ferritic  stainless  steels,  thermodynamic  modeling  was  performed  in  conjunction  with 
experimental  oxidation  testing  and  empirical  kinetic  evaluation.  The  phase  stability  and  oxidation  behav¬ 
ior  of  the  Fe-Cr-0  ternary  system  were  assessed  using  thermodynamic  calculations.  Calculated  ternary 
phase  diagrams  were  validated  against  the  experimental  oxidation  data  of  Fe-20Cr  and  Fe-18Cr  ferritic 
stainless  steels,  GE-13L  and  AL-441HP,  respectively.  Results  indicate  that  the  use  of  accelerated  testing, 
such  as  exposing  the  system  to  higher  temperatures,  can  lead  to  changes  in  phase  equilibria  and  the 
oxidation  kinetics  of  the  alloys.  Through  combined  thermodynamic  assessment  and  controlled  oxidation 
experiments,  the  oxidation  behavior  of  high-chromium  ferritic  stainless  steels  is  presented  and  discussed. 

©  2010  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Solid  oxide  fuel  cells  (SOFCs)  are  solid  state  energy  conversion 
devices  that  operate  by  electrochemically  converting  fuel  and  oxi¬ 
dant  gases  into  electricity  without  combustion.  The  reduction  in 
SOFC  operating  temperature  due  to  recent  development  in  mate¬ 
rials  and  fabrication  techniques  enable  the  use  of  ferritic  stainless 
steels  as  interconnects  for  SOFCs  [1,2].  Ferritic  stainless  steels  offer 
many  advantages  over  traditional  ceramic  interconnects,  such  as 
higher  thermal  and  electronic  conductivity,  reasonable  thermal 
expansion  matching,  relative  ease  of  fabrication  and  decreased 
machining  costs  [3].  However,  the  formation  of  thermally  grown 
chromium  oxide  on  ferritic  stainless  steel  alloys  results  in  an 
increase  in  ohmic  resistance  and  can  lead  to  the  loss  of  cathode 
catalytic  activity  (through  Cr  volatilization)  known  as  cathode  poi¬ 
soning  [4]. 

Oxidation  behavior  of  ferritic  stainless  steels  has  been  exten¬ 
sively  studied.  It  has  been  found  that  oxide  growth  depends  on 
many  factors  including  alloy  composition,  oxygen  partial  pressure, 
temperature,  and  exposure  time  [5-8].  The  effect  of  alloy  com¬ 
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position  on  the  oxidation  behavior  of  engineering  alloys  can  be 
complex,  in  particular  when  dealing  with  minor  alloying  and  impu¬ 
rity  elements.  As  such,  the  current  study  primarily  investigates  the 
effects  of  major  alloying  elements,  specifically  Cr  and  Fe,  on  the 
oxidation  behavior. 

When  exposed  to  air  or  oxygen  at  800-1 000  °C,  ferritic  stain¬ 
less  steels  form  chromium  oxide,  (Fe,Cr,Mn)304-Spinel,  iron 
oxides,  and/or  their  combinations  such  as  iron-chromium  oxides. 
Kubaschewski  and  Hopkins  [9]  proposed  that  the  rate  of  the 
early  oxidation  of  iron-chromium  alloy  is  rapid  due  to  the  for¬ 
mation  of  iron  oxide.  During  the  initial  oxidation  stages,  iron-rich 
Corundum  oxide  (Fe203-Cr203)  forms,  subsequently  converting  to 
chromium-rich  oxide  with  time  [8].  SOFC  performance  is  strongly 
influenced  by  the  electrical  resistance  of  the  stainless  steel  inter¬ 
connect.  Therefore,  increased  cell  resistance  as  a  result  of  the  oxide 
scale  is  a  direct  consequence.  Hence,  it  is  important  to  identify 
phase  stability,  material  responses  under  the  oxidation  conditions, 
and  the  oxidation  kinetics  of  stainless  steel  interconnects  at  oper¬ 
ating  temperatures  (700-850  °C).  Given  anticipated  lifetimes  of 
40,000  h,  accelerated  testing  is  desirable.  However,  accelerated 
testing  at  elevated  temperatures  may  influence  phase  equilib¬ 
ria  and  oxidation  mechanism  and  kinetics  of  the  system  to  the 
point  where  high  temperature  data  is  no  longer  relevant  to  actual 
operating  temperatures.  A  systematic  understanding  of  oxidation 
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behavior  and  oxidation  kinetics  of  ferritic  stainless  steel  intercon¬ 
nects  as  a  function  of  temperature  will  enable  the  prediction  of 
the  long  term  effects  of  chromium  oxide  that  forms  as  a  result  of 
exposure  to  SOFC  cathode  operating  conditions. 

The  aim  of  this  work  is  to  assess  the  oxidation  behavior  of 
ferritic  stainless  steels  using  thermodynamic  modeling  and  exper¬ 
imentally  derived  kinetics.  Two  ferritic  steels  were  selected  for 
this  study:  GE-13L  and  AL-441HP.  GE-13L  is  an  alloy  developed  by 
GE  for  SOFC  interconnect  applications.  AL-441HP  is  a  commercial 
chromia-forming  stainless  steel  available  from  Allegheny  Ludlum 
which  is  a  lower  cost  alternative.  Thermodynamic  modeling  will 
be  presented  in  terms  of  isothermal  sections,  stability  diagrams,  as 
well  as  calculated  phase  fraction  diagrams.  Oxidation  behavior  of 
the  two  steels  will  be  discussed  in  term  of  kinetics  assessment  over 
the  temperature  range  of  800-1 000  °C. 

2.  Methodology 

2.1.  Thermodynamic  modeling 

A  consistent  thermodynamic  description  of  Fe-Cr-0  ternary 
was  obtained  based  on  the  published  thermodynamic  modeling  of 
the  quaternary  Fe-Cr-Ni-0  system  [  1 5  ].  All  of  the  solid  phases  were 
modeled  by  the  Compound  Energy  Formalism  (CEF)  [16],  which 
has  been  widely  used  to  describe  phases  using  two  or  more  sublat¬ 
tices  in  CALPHAD  assessments.  The  molar  Gibbs  energy  of  a  phase 
consists  of  three  terms  as  given  below, 

c£  =  re/c£  +  idc£  +  exc£  (i) 

with  the  first  term  on  the  right  hand  side  of  Eq.  (1)  being  the  Gibbs 
energy  in  the  reference  state,  the  second  term  the  ideal  mixing 
term,  and  the  third  one  the  excess  term,  which  may  be  determined 
through  an  optimization  procedure.  For  this  particular  case,  since 
all  of  the  solid  metallic  phases  and  oxide  phases  undergo  a  mag¬ 
netic  transition,  the  magnetic  contribution  to  the  Gibbs  energy  was 
incorporated  [  1 7,21  ].  The  magnetic  contribution  is  described  by  Eq. 
(2). 

magG  =  RTln(/3  +  1)/(t)  (2) 

where  r  =  T/Tc,  Tc  is  the  Curie  (or  Neel)  temperature  and  is  the 
effective  magnetic  moment  [17,21].  It  is  worth  pointing  out  that 
in  Ref.  [15]  one  of  the  parameters  for  the  Cr-rich  Spinel  phase  was 
found  to  be  in  error,  which  causes  the  Cr304  phase  to  be  artificially 
stable  at  low  temperatures.  A  correction  was  made  to  fix  this  prob¬ 
lem  by  adjust  parameters  for  those  end  members  with  Cr2+  on  the 
first  sublattice. 

In  addition  to  a  reliable  thermodynamic  description,  robust  soft¬ 
ware  is  essential  for  the  successful  application  of  the  CALPHAD 
approach.  In  the  present  study,  Pandat  8.0  [22]  was  used  to  per¬ 
form  the  phase  diagram  calculations.  The  overall  system  pressure 
was  kept  at  1  atm  for  all  calculations. 

2.2.  Alloy  selection  for  oxidation  testing 

GE-13L  is  a  chromia  forming  ferritic  stainless  steel  developed 
by  GE  specifically  for  SOFC  interconnect  applications.  It  is  an 
Fe-20  wt%  Cr  steel  (Table  1)  that  was  optimized  for  machinability, 
strength,  and  oxidation  resistance  at  an  SOFC  operating  tempera¬ 
ture  of  ~800  °C.  Given  that  silicon  and  aluminum  impurities  in  SOFC 
interconnect  alloys  can  lead  to  the  formation  of  highly  resistive 
oxide  scales,  the  processing  of  this  alloy  was  performed  in  vacuum 
to  limit  the  inclusion  of  these  detrimental  elements. 

Stainless  steel  441  (AL-441HP)  is  a  commercial  ferritic  stainless 
steel  with  1 8  wt%  chromium  (Table  1 ).  AL-441  HP  is  readily  available 
and  does  not  require  sophisticated  alloy  manufacturing  methods 


(such  as  vacuum  processing)  making  it  more  cost  effective.  There¬ 
fore,  AL-441  HP  is  an  attractive  option  as  a  lower  cost  option  for 
SOFC  interconnects. 

The  oxidation  study  of  GE-13L  and  AL-441  HP  was  aimed  at 
understanding  the  oxidation  behavior  and  kinetics  of  both  stain¬ 
less  steels  as  a  function  of  temperature  in  the  range  of 800-1000  °C, 
where  800  °C  is  the  target  SOFC  operating  temperature.  However, 
accelerated  testing  through  higher  temperatures  such  as  900  or 
1000°C,  must  be  well  understood  in  order  to  evaluate  the  validity 
of  such  techniques. 

The  oxidation  data  will  be  used  to  identify  any  fundamental  dif¬ 
ferences  in  oxidation  behavior  between  the  two  steels.  The  focus  of 
the  current  study  will  be  to  identify  oxide  phases  formed,  character¬ 
ize  oxidation  behavior,  and  establish  the  parabolic  rate  constants. 
The  resulting  understanding  will  enable  prediction  of  long-term 
degradation  effects  (in  terms  of  oxide  thicknesses)  due  to  oxide 
scale  formation  for  the  two  candidate  materials. 

2.3.  Oxidation  testing 

The  oxidation  behavior  of  ferritic  stainless  steels  GE-13L  and 
AL-441  HP  was  studied  using  cold-rolled  sheet  material.  Samples 
of  25  mm  x  25  mm  x  0.5  mm  or  25  mm  x  12.5  mm  x  0.5  mm  were 
machined  from  each  steel  for  oxidation  tests. 

Individual  samples  were  ultrasonically  cleaned,  weighed,  and 
physically  measured  prior  to  oxidation  testing.  Continuous  isother¬ 
mal  oxidation  tests  were  conducted  at  800,  900,  and  1000°C  in 
laboratory  air  up  to  200  h.  After  testing,  samples  were  removed 
from  the  furnace,  air  cooled  and  weighed  at  room  temperature. 
Weight  change  was  determined  from  the  difference  before  and 
after  the  oxidation  testing  normalized  by  the  sample  surface  area. 
The  weight  gain  of  the  samples  is  indicative  of  oxygen  ingress  and 
formation  of  an  oxide  scale  while  weight  loss  indicates  the  loss  of 
material  through  spallation  and  the  loss  via  volatilization  of  oxides 
such  as  chromia. 

The  oxide  scale  surfaces  were  examined  by  X-ray  diffraction, 
after  which  the  oxidized  samples  were  sectioned  and  polished 
to  characterize  the  oxide  scale  in  a  cross-sectional  orientation. 
The  cross  sections  were  analyzed  by  scanning  electron  microscopy 
(SEM)  equipped  with  energy  dispersive  X-ray  spectroscopy  (EDS). 

2.4.  Kinetic  analysis  methodology 

A  kinetic  model  for  each  testing  temperature  was  determined 
by  Eq.  (3)  [10].  The  kinetic  parameter,  n,  signifies  the  kinetic  model 
for  the  oxidation  behavior  of  materials.  For  the  ideal  diffusion- 
controlled  oxidation  of  pure  metals,  the  oxidation  kinetics  should 
follow  a  parabolic  rate  law  represented  in  Eq.  (4)  [24].  Using  the 
weight  change  and  the  oxide  thickness  data,  the  kinetic  parameter 
(n)  and  parabolic  rates  for  samples  oxidized  at  800, 900,  and  1 000  °C 
were  calculated  and  are  presented  and  discussed  in  the  following 
section. 

Oxidation  kinetic  model  calculation: 


where  Aw  is  the  weight  change,  Area  is  the  surface  area  of  the 
sample,  k  is  a  constant,  t  is  the  time,  and  n  is  the  kinetic  model 
parameter. 

Parabolic  rate  law  calculation: 

W 2  =  kpt  +  C  (4) 

where  W  is  the  weight  change  per  unit  area,  kp  is  a  parabolic  con¬ 
stant,  and  C  is  a  constant. 


Table  1 

Chemistries  of  GE-1 3L  and  AL-441  HP  [  1 8 ]. 
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Alloys 

Weight  (%) 

Fe 

Cr 

Ti 

Nb 

Si 

Mn 

W 

Others 

GE-13L 

Bal. 

20 

4-5 

RE 

AL-441  HP 

Bal. 

18 

0.1-0.5 

1 .0  max 

1 .0  max 

1 .0  max 

C,  P,  S 

Note :  RE  is  reactive  element. 


o 


o 


FE  x%(CR)  CR 

Fig.  1.  Calculated  isothermal  section  of  Fe-Cr-0  at  800  °C. 


Fig.  2.  Calculated  isothermal  section  of  Fe-Cr-0  at  1000  °C. 


3.  Results 

3.1.  Calculated  ternary  isothermal  sections  at  800  and  1 000  °C 

Fig.  1  shows  a  calculated  phase  diagram  of  the  Fe-Cr-0  sys¬ 
tem  at  800  °C.  In  the  present  study,  all  calculations  are  based 
on  atomic  percent  unless  indicated  otherwise.  According  to  this 
isothermal  section,  the  stable  oxides  are  Corundum  [(Fe,Cr)203], 
Fe-rich  Spinel  [(Fe,Cr)304],  and  Halite  [FeO].  The  Corundum  phase 
forms  a  continuous  solid  solution  between  two  binary  com¬ 
pounds:  Fe203  and  Cr203.  In  addition,  the  oxygen  concentration 
in  Corundum  is  constant  at  60at.%.  Given  that  Cr304  is  not  sta¬ 
ble  at  this  temperature,  the  FeCr204-Spinel  phase  extends  from 
the  binary  Fe304  to  FeCr204.  There  is  a  large  two-phase  field 
between  BCC-(Fe)  and  the  Cr-rich  (Fe,Cr)304  Spinel  phase,  Cr203. 
This  indicates  that  for  a  Fe-rich  alloy,  for  example,  Fe-20Cr,  the 
BCC  phase  will  be  stable  at  this  temperature  and  it  should  be  in 
equilibrium  directly  with  Cr203.  The  same  applies  for  the  Cr-rich 
corner. 

Fig.  2  shows  a  calculated  isothermal  section  at  1 000  °C.  The  (0)- 
rich  portion  is  almost  identical  to  the  diagram  at  800  °C,  except 
the  oxygen  solubility  is  slightly  higher  in  the  Spinel  phase.  For  the 
Fe-rich  region,  Cr  concentration  below  ~12  at.%  results  in  FCC-(Fe) 
phase  whereas  BCC-(Fe)  is  a  stable  phase  in  alloys  with  Cr  con¬ 
centration  greater  than  12at.%.  Moreover,  since  the  temperature 
is  above  the  critical  point  of  the  miscibility  gap  between  BCC-(Fe) 
and  BCC-(Cr)  (binary  phase  diagram),  the  two  BCC  phases  become 
undistinguishable.  As  a  result,  the  two-phase  field  of  BCC  +  Cr203  is 
significantly  enlarged.  In  the  current  study,  the  isothermal  section 
for  900  °C  was  also  calculated.  However,  it  is  found  that  the  phase 


equilibrium  relationship  is  very  similar  to  that  of  1000  °C  as  shown 
in  Fig.  2.  Therefore,  it  is  not  shown  here. 

3.2.  Stability  diagram  (predominance,  or  chemical  potential 
diagram ) 

Fig.  3  shows  the  calculated  stability  diagram  at  800  °C.  At  P02 
higher  than  1 0-8 ,  Corundum  is  the  most  stable  phase.  Depending  on 
the  composition  of  the  alloy,  the  Corundum  phase  could  be  either 
Fe-rich  or  Cr-rich.  In  the  intermediate  P02  regime,  several  oxides 


Fig.  3.  Calculated  stability  diagram  (also  called  predominance  diagram,  or  chemical 
potential  diagram)  of  Fe-Cr-0  at  800  °C.  The  blue  dashed  line  indicates  the  Fe-20Cr 
alloy  composition.  X  indicates  atomic  fraction. 
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Fig.  4.  Calculated  stability  diagram  (also  called  predominance  diagram,  or  chemical 
potential  diagram)  of  Fe-Cr-0  at  1000°C.  Note  that  the  blue  dashed  line  indicates 
the  Fe-20Cr  alloy  composition. 


may  present,  including  Spinel,  Halite,  and  Corundum-Cr203 .  At  low 
oxygen  partial  pressures,  BCC  and  the  cr  phase  are  in  equilibrium 
with  Cr-rich  (Cr,Fe)203  Corundum  phase.  For  Fe-20Cr  stainless 
steel,  sequence  of  oxide  formation  is  illustrated  by  the  schematic 
diagram.  It  indicates  that  at  the  oxide/metal  interface  (low  P02 
region)  Cr203  is  anticipated  to  form  in  contact  with  the  metallic 
stainless  steel.  At  the  oxide  surface,  the  phase  stability  diagram  pre¬ 
dicts  Fe203.  In  between  these  two  layers,  several  oxides  may  form, 
i.e.  Fe304,  FeO,  FeCr204-Spinel.  When  compared  to  experimental 
observations,  thermodynamic  calculations  are  in  agreement  with 
oxide  phases  that  are  found  in  the  oxidized  stainless  steels.  Detailed 
experimental  observations  and  thermodynamic  calculations  are 
discussed  later. 

Fig.  4  shows  the  calculated  stability  diagram  at  1000  °C.  It  has  a 
similar  trend  to  that  at  800  °C.  However,  there  are  two  main  differ¬ 
ences.  First,  at  the  oxide/metal  interface,  FCC-(Fe)  may  form  instead 
of  the  a  phase  as  in  the  800  °C  case.  The  implication  of  the  existence 
of  the  FCC  phase  is  that  the  diffusion  kinetics  may  be  altered.  Sec¬ 
ond,  each  phase  equilibrium  window  shifts  toward  higher  P02.  As 
in  the  figure  for  800  °C,  the  oxide  formation  sequence  is  represented 
by  the  schematic  diagram. 

3.3.  Oxidation  behavior  of  ferritic  steels  as  a  function  of 
temperatures 

The  oxidation  behavior  of  GE-13L  and  AL-441HP  were  mea¬ 
sured  by  changes  in  weight  and  oxide  thickness  as  a  function  of 
temperature  and  time.  Weight  changes  of  the  steels  were  normal¬ 
ized  to  the  original  area  and  are  summarized  in  Fig.  5.  All  samples 
gained  weight  throughout  the  test  indicating  oxygen  ingress  and 
the  formation  of  oxides.  For  GE-13L,  iron-chromium  oxides  are  the 
primary  oxide.  For  AL-441HP,  weight  gain  was  not  only  a  result 
of  faster  oxide  growth  rate  the  formation  of  iron-chromium  and 
iron-chromium-manganese  Spinel  oxides,  but  also  a  result  of  the 
internal  oxidation  of  titanium. 

Initial  weight  gain  of  both  GE-13L  and  AL-441HP  was  high 
followed  by  lower  weight  change  rate  as  protective  oxide  scale 
formed.  However,  rate  of  change  for  AL-441HP  was  higher  than 
that  of  GE-13L.  In  other  word,  AL-441HP  continued  to  gain  weight 
at  a  faster  rate  than  GE-13L  did.  Moreover,  the  total  weight  gain  of 
AL-441  HP  was  in  an  order  of  magnitude  greater  than  that  of  GE-1 3L 
(depending  on  the  temperature).  Higher  weight  gain  of  AL-441  HP 
was  a  result  of  both  faster  oxidation  kinetics  (as  shown  by  thicker 
oxide  scale)  and  internal  oxidation  of  Ti02. 

In  order  to  determine  the  phases  formed  during  oxidation,  the 
surfaces  of  the  oxidized  samples  were  characterized  by  X-Ray 
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Fig.  5.  Comparison  of  weight  change  as  a  function  of  time  and  temperature  for  GE- 
13L  (circles)  and  AL-441  HP  (triangles).  The  parabolic  fits  are  shown  in  the  graph 
with  solid  lines  (for  GE13L)  and  dashed  lines  (for  AL-441  HP).  R 2  values  indicating 
the  fit  of  data  are  shown  along  with  each  parabolic  fit. 


diffraction.  Three  phases  were  detected  on  oxidized  GE-1 3L:  Corun¬ 
dum  chromium-iron  oxide  [(Cr,Fe)203],  Spinel  chromium-iron 
oxide  [(Cr,Fe)304],  and  iron  [Fe].  The  detection  of  base  metal,  Fe, 
indicates  that  the  oxide  scale  was  thin  (<25  pm),  allowing  the  X- 
ray  interaction  volume  to  penetrate  through  the  oxide  layer  and 
sample  the  bulk  alloy.  For  AL-441  HP,  XRD  of  the  oxidized  samples 
revealed  all  three  phases  detected  in  oxidized  GE-13L.  However,  an 
additional  oxide  of  rutile  (Ti02)  was  also  detected.  Another  main 
difference  between  GE-1 3L  and  AL-441  HP  was  the  peak  height  ratio 
of  Corundum/Spinel  as  shown  in  Table  2.  The  Corundum/Spinel 
peak  ratio  in  both  alloys  decreased  with  time  and  temperature  sug¬ 
gesting  higher  volume  fraction  of  Spinel  (Cr,Fe,Mn)304  in  the  oxide 
scale.  When  compared  between  alloys,  GE-13L  showed  higher 
Corundum/Spinel  peak  height  ratio  which  indicated  higher  volume 
fraction  of  Corundum  (Cr,Fe)203  in  GE-13L 

Experimental  results  indicate  the  solubility  of  iron  and 
chromium  in  both  Corundum-Cr203  (hematite-type  crystal  struc¬ 
ture)  and  Spinel-(Fe,Cr,Mn)304  phase  (magnetite-type  crystal 
structure),  in  good  agreement  with  thermodynamic  calculations. 


Table  2 

Relative  peak  height  ratio  of  Corundum  (Cr,Fe)203  and  Spinel  (Fe,Cr,Mn)304  in  GE- 
13L  and  AL-441HP  after  oxidation  at  800,  900,  and  1000  °C.  Phases  were  identified 
by  X-ray  diffraction  (XRD). 


Temperature 

Testing  time 

Alloys 

Corundum/Spinel 
peak  height  ratio 

800  °C 

25  h 

GE-13L 

1.0 

AL-441  HP 

1.0 

100  h 

GE-13L 

1.0 

AL-441  HP 

0.6 

900  °C 

6h 

GE-13L 

1.0 

AL-441  HP 

0.8 

25  h 

GE-13L 

1.0 

AL-441  HP 

0.7 

100  h 

GE-13L 

0.9 

AL-441  HP 

0.3 

1000 °C 

6h 

GE-13L 

0.9 

AL-441  HP 

0.5 

25  h 

GE-13L 

0.9 

AL-441  HP 

0.1 
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Fig.  6.  Comparison  of  oxide  thickness  as  a  function  of  temperature  and  alloys,  GE-13L  vs.  AL-441HP. 


The  main  difference  between  the  experimental  observations  and 
the  thermodynamic  calculations  relates  to  the  outer  oxide  layer. 
The  stability  diagrams  at  both  800  °C  (Fig.  3)  and  1000°C  (Fig.  4) 
indicate  that,  as  a  function  of  oxygen  partial  pressure,  Fe203  oxide 
should  form  as  the  outer  layer,  followed  by  Spinel-(Fe,Cr,Mn)304 
while  Cr203  oxide  forms  as  a  sublayer  at  the  oxide/metal  inter¬ 
face  for  both  GE-13L  and  AL-441HP.  Flowever,  experimental  data 
showed  only  iron-chromium  oxide  scales  of  Corundum  Cr203- 
and  Spinel  (Fe,Cr,Mn)203-type.  The  discrepancy  could  be  due  to 
the  kinetics  of  the  oxide  formation  involving  the  diffusion  of  each 
species,  the  time  required  for  reactions  to  equilibrate,  the  differ¬ 
ence  of  partial  pressure  of  oxygen,  or  the  loss  of  the  outer  oxide 
layer  during  sample  mount  preparation.  Despite  the  discrepancy, 
the  trend  of  oxide  formation  (iron-chromium  oxide  scale)  was  self- 
consistent. 

After  surface  analysis  by  XRD,  further  characterization  was 
performed  by  examining  cross  sectional  samples  using  an  SEM 
equipped  with  EDS.  Back  scattered  electron  images  of  selected 
cross-sectioned  oxidized  samples  are  shown  in  Fig.  6.  Oxide  thick¬ 
nesses  were  measured  directly  from  cross  sections  of  oxidized 
alloys  and  are  summarized  in  Table  3.  Oxide  thickness  data  corre¬ 
lates  well  with  weight  change  data  and  no  spallation  was  observed 
at  all  testing  temperatures  through  test  durations. 

Back-scattered  electron  images  of  cross  section  of  oxidized  GE- 
13L  showed  a  thin  and  dense  iron-chromium  oxide  layer.  Scratches 
evident  in  the  cross  section  images  are  due  to  metallographic 
effects  (sectioning  and  polishing  of  the  specimens).  Energy  disper¬ 
sive  spectroscopy  (EDS)  was  used  to  detect  elements  present  within 
different  layers.  The  EDS  analyses  revealed  that  the  oxide  layer  con¬ 


sisted  primarily  of  chromium  and  oxygen.  A  small  peak  of  Fe  was 
also  detected  potentially  indicating  the  solubility  of  Fe  in  chromium 
oxide  or  possibly  Fe-Cr  Spinel.  Both  phases  are  thermodynamically 
stable  as  indicated  earlier  in  thermodynamic  model.  Moreover,  due 
to  the  thin  oxide  scale,  the  interaction  volume  of  the  electron  beam 
could  sample  Fe  from  base  alloy. 

For  AL-441HP,  the  oxide  scale  of  the  800  °C  sample  was  dense 
and  thin  (1.1  [xm  after  lOOh  at  800  °C).  Oxide  scales  at  900  and 
1 000  °C  were  dense  similar  to  that  of  GE-1 3L.  EDS  analysis  of  oxide 
scale  indicated  that  the  oxide  layer  consisted  mostly  of  chromium 
and  oxygen  with  a  small  amount  of  iron.  Internal  oxidation  primar¬ 
ily  consisted  of  Ti02  was  observed  in  900  °C  and  1000°C  samples. 
The  depth  of  internal  oxidation  increased  with  increasing  temper¬ 
ature.  AL-441HP  contains  between  0.1  and  0.6  wt%  Ti  and  internal 
oxidation  of  a  stainless  steel  containing  titanium  was  also  observed 
by  Fluczkowski  et  al.  [19].  Internal  oxidation  of  Ti02  contributed, 
in  part,  to  higher  weight  gain  observed  in  AL-441FIP  (compared  to 
GE-13L). 

4.  Discussion 

4. 1 .  Kinetic  analysis  at  800 0  C,  900 0  C,  and  1 000 0  C 

The  rate  of  weight  change  as  a  function  of  temperature  is  signif¬ 
icantly  greater  for  AL-441HP  when  compared  to  GE-13L  at  800  °C, 
900  °C,  and  1000°C.  Both  the  internal  oxidation  and  the  lower 
chromium  content  in  AL-441HP  likely  contribute  to  the  higher 
weight  gain  as  well  as  the  higher  rate  of  change  due  to  higher  vol¬ 
ume  fraction  of  Spinel-(Cr,Fe)304  (Table  2)  non-protective  iron-rich 


Table  3 

Oxide  thickness  of  GE-13L  and  AL-441HP  oxidized  at  800,  900,  and  1000  °C. 

Testing  time  Oxide  thickness  (|jim) 

800  °C  900 °C  1000°C 


GE-13L  AL-441HP  GE-13L  AL-441HP  GE-13L  AL-441HP 

6h  Nil  Nil  0.7  1.4  1.2  3.9 

25  h  0.2  0.8  1.1  1.8  3.1  8.3 

lOOh  0.5  1.1  1.7  3.5  8.6  13.1 


1980 
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Fig.  7.  The  oxidation  rate  constant  as  a  function  of  inverse  temperature  of  GE-13L 
and  AL-441HP  oxidized  at  800-1000  °C. 

oxide  scale.  No  internal  oxidation  was  observed  in  GE-13L  at  any 
test  temperature. 

The  oxidation  behaviors  of  GE-13L  and  AL-441HP  follow  the 
parabolic  rate  law.  The  oxidation  kinetics  were  assessed  based  on 
weight  change  data  and  are  shown  in  the  Arrhenius  plot  in  Fig.  7. 
Parabolic  constants  of  both  GE-13L  and  AL-441HP  increase  with 
temperature  indicating  enhanced  kinetics  at  higher  temperature. 

The  parabolic  rate  constants  of  AL-441  HP  are  an  order  of  magni¬ 
tude  greater  than  that  of  GE-13L  at  all  temperatures.  As  such,  oxide 
scales  on  AL-441  HP  are  consistently  thicker  than  that  on  GE-13L  at 
similar  times  and  temperatures  (Fig.  6).  The  faster  kinetics  of  AL- 
441  HP  oxidation  is  likely  a  combined  result  of  internal  oxidation 
and  lower  chromium  content.  The  effect  of  the  chromium  content 
will  be  discussed  in  the  next  section  (Section  4.2)  by  combining  the 
experimental  results  with  thermodynamic  calculations. 

The  oxidation  behavior  of  AL-441  HP  is  more  complex  than 
GE-13L  due  to  impurities/alloying  elements  present  in  AL-441  HP. 
Although  a  chromium  oxide  former,  AL-441  HP  contains  titanium 
as  well  as  silicon,  both  of  which  are  oxide  formers.  Internal  oxida¬ 
tion  of  AL-441  HP  was  evident  in  this  oxidation  study  as  indicated 
in  XRD  and  cross  section  results  (Fig.  6).  The  formation  of  silica 
was  not  detected  by  the  methods  used  in  the  examination  of  oxi¬ 
dized  AL-441  HP.  However,  the  test  duration  may  not  be  sufficient 
for  silica  scale  to  form.  Jablonski  et  al.  [20]  reported  the  formation 
of  silica  in  samples  oxidized  at  800  °C  for  2000  h  by  TEM.  Therefore, 
further  testing  and  more  detailed  characterization  should  be  done 
to  ensure  the  absence  of  continuous  silica  scale  under  long  term 
operation. 

4.2.  Effects  of  chromium  concentrations 

Chromium  content  is  expected  to  be  a  primary  factor  in  the  for¬ 
mation  of  chromium  oxide  scale  in  ferritic  stainless  steels.  In  this 
study,  AL-441  HP  with  18  wt%  Cr  is  compared  to  GE-13L  which  has 
20  wt%  Cr.  In  order  to  investigate  effects  of  chromium  content  in  this 
range,  thermodynamic  calculations  were  performed  and  compared 
with  experimental  data. 

Fig.  8(a)  and  (b)  shows  the  calculated  phase  fractions  for  the 
Corundum  and  Spinel  phases  as  a  function  of  P02  in  log  scale  for 
Fe-20Cr  and  Fe-1 8Cr  alloys  at  900  °C  and  1 000  °C,  respectively.  For 
both  alloys,  with  increasing  Cr  content,  the  relative  fraction  of  Cr- 


Fig.  8.  Calculated  the  fraction  of  Corundum  and  Spinel  phases  as  a  function  of  P02 
in  two  alloys:  Fe-20Cr  and  Fe-18Cr:  (a)  900°C  and  (b)  1000°C. 


rich  Spinel  ((Cr,Fe)304)  and  Cr-rich  Corundum  (Cr203)  increases. 
However,  the  volume  fraction  of  the  Fe-rich  Corundum  (Fe203 )  and 
the  Fe-rich  Spinel  ((Fe,Cr)304)  is  independent  of  the  Cr  content. 
Furthermore,  as  shown  in  Table  4,  the  calculations  also  indicate 
that  for  both  temperatures  (900  and  1000°C),  the  Fe-20Cr  tends 
to  have  a  slightly  higher  ratio  of  Corundum  to  Spinel  content  than 
the  Fe-18Cr.  The  trend  of  higher  Corundum/Spinel  ratio  was  also 
observed  in  the  oxidation  experiment  (Table  2). 

Thermodynamic  calculations  suggest  that  higher  chromium 
content  results  in  an  oxide  with  a  higher  fraction  of  Corundum 
((Cr,Fe)203)  compared  to  the  Spinel  ((Fe,Cr)304)  phase.  This  result 
is  consistent  with  the  experimental  results  (Table  2).  At  900  °C, 
the  oxide  scale  of  GE-13L  (20wt%  Cr)  appears  to  have  a  higher 
weight  fraction  of  Corundum  phase  (Table  2)  than  AL-441  HP.  At  all 
testing  conditions  (temperature  &  time),  the  relative  ratio  of  Corun¬ 
dum/Spinel  phase  in  GE-13L  is  higher  than  in  AL-441  HP  (Table  2). 
It  should  be  pointed  out  that  Mn  which  is  present  in  AL-441  HP  but 


Table  4 

Calculated  ratio  between  the  fraction  of  Corundum  and  that  of  Spinel  phase  at  a 
specific  P02  where  these  two  phases  are  in  equilibrium. 


Condition 

Phase  amount  ratio:  Corundum/Spinel 

Fe-20Cr 

Fe-18Cr 

900  °C  (log  P02  =  -20.58) 

0.771 

0.767 

1000  °C  (log  P02  =  -18.37) 

0.774 

0.770 
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Fig.  9.  Calculated  Cr-content  in  the  Corundum  phase  as  a  function  ofP02  in  log  scale 
for  Fe-20Cr  and  Fe-18Cr  alloys,  (a)  900  °C  and  (b)  1000  °C. 


not  in  GE-13L  is  a  Spinel-phase  stabilizer.  Hence,  the  presence  of 
Mn  could  decrease  the  Corundum/Spinel  ratio  in  AL-441HP. 

Since  Cr  and  Fe  are  miscible  in  both  the  Corundum  and  Spinel 
phases,  the  difference  in  chromium  concentration  may  affect  not 
only  the  fraction  of  oxides  formed,  but  also  the  ratio  of  Cr/Fe  in 
each  phase.  Fig.  9  shows  the  calculated  Cr-content  in  the  Corun¬ 
dum  phase  as  a  function  of  PO2  in  log  scale  for  Fe-20Cr  and  Fe-1 8Cr 
alloys  at  (a)  900  °C  and  (b)  1 000  °C.  For  both  temperatures,  the  cal¬ 
culations  show  that  in  the  Fe-rich  Corundum  phase,  the  Cr-content 
appears  to  follow  the  same  trend  as  the  overall  Cr  content  in  the 
alloy.  In  other  words,  the  Fe-rich  Corundum  phase  (Fe,Cr)203  in 
Fe-20Cr  has  slightly  higher  Cr-content  than  that  of  Fe-18Cr.  How¬ 
ever,  the  Cr-content  in  the  Cr-rich  Corundum  phase  (Cr,Fe)203 
seems  to  be  insensitive  to  the  overall  Cr  concentration  in  the  alloy. 

As  previously  shown  in  Fig.  6  when  the  oxidation  behaviors  of 
the  two  steels  are  compared  under  identical  conditions,  it  is  clear 
that  the  oxide  scale  of  GE-1 3L  is  thinner  than  that  of  AL-441  HP,  sug¬ 
gesting  a  more  protective  oxide  scale.  Kubaschewski  and  Hopkins 
[9]  proposed  that  in  chromia-forming  steels,  oxides  with  higher 
chromium  concentration  in  the  Corundum  and  Spinel  phases  are 
more  desirable  due  to  increased  protectiveness.  However,  this 
interpretation  takes  into  account  only  the  chromium  content  and 
not  the  effects  of  other  alloying  elements.  In  AL-441  HP,  other  alloy¬ 
ing  elements  such  as  Ti  form  oxides  that  affect  the  weight  change  of 
the  oxidized  samples.  Gertsricken  and  Dekhtyar  [26]  investigated 
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Fig.  10.  Comparison  of  oxide  thicknesses  between  calculated  (est,  solid  and  dash 
lines)  and  experimental  results  (exp,  circle  and  triangle  symbols)  of  GE-13L  and 
AL-441  HP  oxidized  at  800  °C. 


the  effects  of  alloying  elements  on  the  diffusion  rate  of  chromium 
in  iron  at  950-1 050  °C.  They  indicated  that  Ti,  Si,  and  Nb  decrease 
chromium  diffusion  rate  in  iron  by  a  half  to  one  order  of  magnitude 
[26].  The  Ti,  Si,  and  Nb  presented  in  AL-441  HP  (but  not  in  GE-13L) 
may  lower  the  kinetic  rate  of  chromium  oxide  formation,  and  thus 
have  a  detrimental  effect  on  the  formation  of  a  protective  oxide 
scale.  Moreover,  the  presence  of  several  alloying  elements  not  only 
individually  influence  the  diffusion  rates  of  Cr  and  Fe,  but  may  also 
have  interaction  effects  on  the  overall  oxidation  behavior  of  the 
stainless  steel. 

For  steels  with  higher  chromium  contents,  thermodynamic  cal¬ 
culations  suggest  increased  Corundum/Spinel  ratios,  which  is  in 
good  agreement  with  experimental  results.  Alloys  with  higher  Cr 
content  also  form  oxides  with  higher  Cr/Fe  ratio  in  each  phase.  GE- 
13L  (20  wt%  Cr)  shows  better  oxidation  resistance  over  AL-441  HP 
(18  wt%  Cr)  from  experimental  results. 

4.3.  Oxide  thickness  analysis 

A  kinetic  analysis  of  GE-13L  and  AL-441  HP  reveals  that  the 
oxidation  of  both  alloys  follows  the  parabolic  rate  law.  The  differ¬ 
ence  between  the  parabolic  constants  of  GE-13L  and  AL-441  HP  is 
substantial  with  parabolic  constants  of  AL-441  HP  an  order  of  mag¬ 
nitude  greater  than  that  of  GE-13L  (Fig.  7).  This  difference  is  also 
evident  in  higher  weight  gain  and  thicker  oxides  present  for  AL- 
441  HP  oxidized  at  800-1 000  °C.  The  higher  parabolic  constant  in 
AL-441  HP  was  determined  from  weight  change  data  (Fig.  5).  High 
weight  change  was  a  result  of  faster  kinetics  as  well  as  the  forma¬ 
tion  of  Ti02  by  internal  oxidation  in  AL-441  HP  (Fig.  6).  AL-441  HP 
could  contain  up  to  0.5  wt%  of  Ti  with  a  typical  concentration  of 
0.2  wt%  [18].  Given  the  amount  of  Ti  in  the  alloy  and  the  internal 
oxidation  depth,  it  is  estimated  that  the  formation  of  Ti02  could 
contribute  1 .3  mg  cm-2  to  weight  gain  at  1 000  °C  for  6  h  (assuming 
all  Ti  oxidized). 

According  to  Wagner’s  theory  of  oxidation  [29,30],  oxide  thick¬ 
ness  of  scale  dominated  by  diffusion  follows  parabolic  rate  law  as 
shown  in  Eq.  (5).  Fig.  10  reports  calculated  results  at  800  °C  along 
with  experimentally  measured  oxide  thickness  data  for  compari¬ 
son.  The  calculated  and  experimental  data  for  GE-13L  are  within 
±5%.  For  AL-441  HP,  the  discrepancy  is  in  the  order  of  +/-  from  0 


1982 


V.S.  Dheeradhada  et  al.  /  Journal  of  Power  Sources  196(2011)  1975-1982 


to  30%  and  may  result  from  the  transition  in  oxidation  behavior 
from  logarithmic  to  parabolic  at  the  early  stages  of  oxidation.  In 
addition,  the  error  associated  with  measuring  thin  oxide  scales  and 
small  weight  change  is  compounded  by  the  assumptions  implicit 
in  the  calculations.  As  exposure  time  increases,  the  discrepancy 
decreases.  For  AL-441HP,  the  discrepancy  is  also  affected  by  inter¬ 
nal  oxidation,  by  a  variation  in  oxide  thickness  and  weight  change 
measurements. 

Parabolic  law  [29]: 


where  §  is  the  oxide  thickness,  kp  is  a  parabolic  rate  constant  in 
weight,  t  is  the  time,  x  is  the  weight  fraction  of  oxygen  in  the  oxide 
(for  example,  Cr203,  x  is  48/152),  and  p  is  the  density. 


5.  Summary 

This  study  contributes  to  the  understanding  of  the  oxidation 
behavior  of  ferritic  stainless  steels  for  SOFC  applications.  A  com¬ 
bination  of  thermodynamic  modeling  and  oxidation  experiments 
provides  a  unique  perspective  on  the  understanding  of  the  oxida¬ 
tion  behavior  of  ferritic  stainless  steels. 

A  self-consistent  thermodynamic  description  (database)  of 
Fe-Cr-0  ternary  system  focusing  on  the  SOFC  operation  tem¬ 
peratures  800-1 000  °C  was  obtained  and  several  observations 
applicable  to  ferritic  stainless  steels  are  found.  First,  on  the  alloy 
side,  FCC-Fe  may  form  at  900  and  1000°C,  with  potential  implica¬ 
tions  on  the  kinetics  of  oxide  scale  growth.  Second,  phase  equilibria 
at  high  temperature  (1000  °C)  shift  toward  higher  P02  region.  Fur¬ 
thermore,  the  amount  and  composition  of  the  oxide  phases  as  a 
function  of  P02  were  calculated  based  on  the  assumption  that  phase 
equilibrium  prevails  in  the  oxidation  process  (which  may  not  be 
the  case).  Results  from  the  thermodynamic  models  were  compared 
with  experimental  data  and  found  to  be  in  good  agreement. 

It  was  also  found  that  the  alloy  Cr  content  will  not  have  sig¬ 
nificant  effect  on  the  oxide  formation  but  does  have  an  appreciable 
influence  on  the  relative  fraction  of  some  of  the  oxide  phases.  Com¬ 
bining  isothermal  sections  and  the  stability  diagrams,  a  possible 
phase  forming  sequence  was  deduced  for  the  oxidation  process  of 
the  SOFC  cathode  interconnect  stainless  steel  alloy.  The  effect  of 
chromium  content  is  evaluated  in  terms  of  the  volume  fraction  of 
two  different  oxides:  Corundum-(Fe,Cr)203  and  Spinel-(Fe,Cr)304. 
It  is  found  that  steel  with  higher  chromium  content  (GE-13L) 
showed  higher  Corundum/Spinel  fraction  ratio.  Furthermore,  alloy 
with  higher  Cr  content  appears  to  result  in  Fe-rich  Corundum  phase 
(Fe,Cr)203  with  higher  Cr  content.  These  results  lead  to  a  more 
protective  oxide  scale  in  alloy  with  higher  Cr  content. 

The  kinetic  evaluation  of  the  thermally  grown  oxide  on  GE-13L 
(20wt%  Cr)  and  AL-441FIP  (18wt%  Cr)  was  performed  and  found 
to  follow  the  parabolic  rate  law  over  the  temperatures  of  interest 
(800-1000  °C).  The  parabolic  constants  of  AL-441  HP  were  an  order 
of  magnitude  higher  than  that  of  GE-13L  at  all  temperature  indi¬ 
cating  faster  kinetics  in  AL-441  HP.  Further  investigation  of  oxide 
morphology  revealed  the  formation  of  thicker  oxide  scale  in  AL- 


441  HP  than  that  in  GE-13L.  Moreover,  internal  oxidation  of  Ti02  in 
AL-441  HP  was  observed  whereas  no  evidence  of  internal  oxidation 
was  observed  in  GE-13L. 

An  oxide  thickness  prediction  model  was  constructed  for  GE- 
13L  and  AL-441  HP  in  the  temperature  ranges  800-1 000  °C.  The 
predicted  results  are  in  good  agreement  with  experimental  data. 

Future  work  includes  a  more  detailed  study  on  the  implications 
of  accelerated  oxidation  testing  though  increased  temperature  to 
better  understand  the  potential  impact  of  the  FCC  phase. 
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